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摘 要 : 针对 低 阶 协调 有 限 元 对 Qi 一 Po, P, 一 PP， 对 二 维 定常 不 可 压缩 Navier-Stokes 方 程 ， 提 出 了 建 
立 在 局 部 压力 投影 上 的 一 类 稳定 化 有 限 元 方法 。 与 其 它 的 稳定 化 方法 相 比 ， 稳 定 项 不 需要 介入 稳 
定 化 参数 ， 不 用 进行 高 阶 导数 的 运算 ， 或 者 边界 积分 ， 稳 定 在 局 部 单元 上 进行 ， 且 容易 编程 实 
现 。 针 对 稳定 化 有 限 元 逼近 解 ， 证 明了 最 优 的 误差 估计 。 数 值 实 验 表 明 ， 该 方法 有 很 好 的 稳定 性 


和 收敛 性 。 
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1 引言 
对 于 粘性 不 可 压缩 流动 的 Navier-Stokes 方 程 ， 其 混合 有 限 元 方法 的 研究 一 直 是 一 个 热点 问 


题 。 为 了 获得 收敛 稳定 的 有 限 元 逼近 ， 通 常 的 速度 /压力 空间 组 合 要 求 满足 Babugka-Brezzi 稳 
定性 条 件 (下 面 简称 B-B 条件)。 众 所 周知 ， 低 阶 协调 元 记 — Po (速度 线性 ， 压 力 常 数 ) 三 角 元 
MQ 一 Po (速度 分 片 双 线性 ， 压 力 分 片 常数 ) 四 边 形 元 是 不 稳定 的 。 虽 然 它 们 不 满足 B-B 条 
件 ， 但 是 ， 它 们 有 着 简单 的 、 正 则 的 数据 结构 ， 而 且 ， 从 计算 机 实现 的 角度 来 看 ， 所 解决 的 三 
维 问题 代数 方程 组 规模 小 ， 带 宽窄 。 因 此 ， 低 阶 协调 元 空间 在 工程 实际 应 用 上 仍 受 到 人 们 的 关 
注 。 

为 了 去 掉 B-B 条 件 的 限制 ， 针 对 低 阶 协调 有 限 元 对 ， 有 限 元 专家 提出 了 很 多 的 稳定 化 方 
法 [9]。 但 是 ， 有 些 方法 的 结果 依赖 于 稳定 化 参数 的 选取 ， 有 些 方 法 还 涉及 到 压力 的 高 阶 导 
数 ， 边 界 积分 和 网 格 的 相互 租 套 等 ， 比 如 文献 [5,6] 针 对 Q, — 马 元 提出 了 局 部 /全 局 的 压力 跳 
跃 形式 ， 此 种 方法 通过 添加 压力 在 元 素 边界 上 的 跳跃 项 来 克服 B-B 条 件 的 限制 ， 但 是 ， 对 于 给 
定 的 网 格 剖 分 ， 很 难 选取 最 优 的 稳定 化 参数 。 

最 近 ，Bochev 等 人 针对 低 阶 P, /P,, Q1/Q1, P, / Po, Qi/Po Tt, 给 出 了 Stokes 方 程 的 压力 
投影 稳定 化 方法 的 理论 分 析 ， 这 种 方法 不 同 于 现 有 的 稳定 化 方法 ， 稳 定 项 不 需要 介入 稳定 化 参 
数 ， 不 用 进行 高 阶 导 数 的 运算 ， 或 者 边界 积分 ， 稳 定 在 局 部 单元 上 进行 ， 且 容易 编程 实现 。 本 
文 主要 针对 Q 一 Po, Pi 一 Po 元 ， 将 这 种 压力 投影 稳定 化 方法 推广 到 定常 Navier-Stokes 方程 ; 
推导 出 最 优 的 误差 估计 ; 数值 实验 表明 ， 该 算法 有 很 好 的 稳定 性 和 收敛 性 。 
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2 定常 Navier-Stokes 方程 的 基本 理论 


设 Q 是 R?2 中 的 有 界 区 域 ， 有 Lipschitz 连续 边界 T。 粘 性 不 可 压缩 流体 的 的 运动 规律 可 由 
如 下 的 Navier-Stokes 方程 


—vAu + (u: V)ju+Vp= f, VrzeQ, (1) 
div =0, Vren, ujr=0 (2) 


进行 描述 ， 其 中 必 = (wi, wuz) 表示 速度 向 量 ，p 是 压力 ， 是 外 力 密度 ，v > 0 是 粘性 系数 。 
为 了 引进 Navier-Stokes 方 程 的 变 分 形式 ， 引 进 下 列 Sobolev 空间 


X = HQP, Y = L2(Q)2, D(A) = B2(Q)2 n X, 


M = (9) = [a e O); f a(z)az = 0 ), 
Q 


六 上 装备 有 通常 的 内 积 和 范 数 
((u, v)) = (gradu, gradv), |u| = ((u,u))”?, 


其 中 (小 ||: Ip Has L2(Q), LO) 上 的 内 积 和 范 数 , | .| 是 Sobolev Ži Hi(Q) a Hi(Q)2 (i = 
1,2) 中 的 范 数 。 由 于 ASQ) 中 的 半 范 和 全 范 是 等 价 的 ， 下 面 文章 中 的 半 范 与 全 范 都 用 | h 
# X x X fl X xM 上 分 别 定义 连续 的 双 线 形 形 式 


a(u,v)=v((u,v)), Vu €X, b(u,q) = (divu,q), vveX, geM 
M(X x M)x(XxM) 上 广义 的 双 线 性 形式 


B((u,p); (u, q)) = a(u,u) — b(v, p) + b(u,q), Y (u, p), (v, q) € (X, M), 


引进 三 线性 形式 
b(u,v,w) = | (u-V)u-wdz, Vu,v,w € X, 
J 
并 给 出 以 下 列 估计 式 
b(u,v,w) = —b(u,w,v), VueVv, xw,u € X, (3) 
bluo, w) < Zeollu||Š lul? (lollo roll? + olè loloh), Yuv wex, (4) 


[6(u, v, w)| + [b(v, u, w)| + [b(w, u, v)| < cilullillvllzllvllo, Yiu € X, ve D(A), u € Y, (5) 
庞 加 莱 不 等 式 成 立 
llvllo < Yollvih, (6) 
这 里 ci = 0,1, p 是 依赖 于 的 常数 。 下 面 文章 中 出 现 的 常数 都 依赖 于 (v, Q, 了)， 不 同 的 
地 方 代表 的 值 不 同 。 
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对 任意 给 定 的 fe Y, (1)-(2) 的 变 分 形式 是 : 求解 (w,p) e X x M， 使 得 
B((u,p); (v,9)) +b(u, u,v) = (f,0), V (u,q) € (X, M). 
E2158 假设 v > 0 和 feY 满足 下 面 的 唯一 性 条 件 
Co 


2 
YO 
> 0, 
v2 If llo 


则 变 分 问题 (7) 有 唯一 的 解 (wp) e (D(A) n X, H1(Q) n M) EWE 


pues 


lula < 2 fllos lulls + Dl < llo. 


3 ”稳定 化 有 限 元 算法 


对 任意 的 h > 0， 令 7 为 人 Q 将 划分 成 闭 的 三 角形 或 四 边 形 KK 的 一 个 剖 分 ， 设 其 满足 通 


正则 性 条 件 ， 即 存在 常数 o > 1 和 0 << w < 1， 使 得 


hk < opk, [cos@k| <w, i=1,2,3,4, VK etm. 


(7) 


(8) 


(9) 


这 里 hx 为 单元 天 的 直径 ，px 为 单元 天 MADAMA, 6 为 四 边 形 单元 剖 分 情形 下 的 天 内 
角 。 定 义 网 格 参数 = max hx, 定义 全 部 单元 之 间 的 内 边界 所 组 成 的 集合 为 Th,， 并 赋予 如 下 


的 范 数 l 
lulle, = uds ) . 
s= (Xf as) 
定义 本 文中 的 有 限 元 子 空间 


P\(K KRZAK 
(K) = (K), 如 果 五 是 fe 
Q (K), WRK ROWE. 


定义 连续 的 分 片 双 线性 速度 子 空间 
Xn={vEX:ivilk ER(K), i=1,2, VK eT,), 
和 分 片 常数 压力 子 空 间 
Mn={q€M:q|x € (K), V K € rh ), 


RU X, 和 ,分别 是 速度 和 压力 函数 空间 XX 和 YY 的 有 限 元 子 空间 。 
于 是 下 面 经 典 的 逼近 结果 成 立 四 


lv — Inullo + hš lv — Isu||r, + h||u — Invlli < Ap2llola，YvweD(4)， 
lla — Jaqllo < «hllalli, V q e H'(Q)n M. 


JET In : D(A) 一 Xn RBA f, Jn: H'(Q) O M > Mp Æ D ERE. 


(10) 


(11) 
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为 了 弥补 低 阶 协调 有 限 元 对 inf-sup 条 件 的 不 足 ， 定 义 如 下 双 线 性 形式 的 稳定 化 项 
G(p,q) = (p — Ip, q — Iq), (12) 


HPI: L2(Q) 一 R, 是 局 部 压力 投影 。 
定常 Navier-Stokes 方 程 建立 在 局 部 压力 投影 上 的 稳定 化 有 限 元 方法 如 下 


Bn ((un, Pa); (v, 9) + b(un, Un, v) 一 (f,v), (13) 
其 中 双 线性 形式 
Bn ((uh, ph); (v, q)) 一 a(un, v) = d(v, pn) =F d(un, q) + G(Ph, q), V (v, q) € (Xh, Mhn). (14) 


下 面 给 出 Stokes 问题 建立 在 局 部 压力 投影 上 的 稳定 化 有 限 元 方法 的 稳定 性 。 
E310] $ (Xn, Mp) 是 最 低 阶 的 协调 有 限 元 对 。 假 设 


jpllo < xlplo，YpeZQ)， lp~Ipllo < «hllplhh, Yp EHQ), (15) 
成 立 ， 则 存在 与 无 关 的 常数 6 > 0, WE 

IBn((u,p); (o,q))| < (llull + iplo ) (lvli + llallo), Y (wu,p),(v,9) E€ (X, M), (16) 

B(llualla + llpallo ) 


Bhx((un, Ph); (Uh, 
< smp Enp) Dl (an, ph), (onah) € (Xr Mn). (17) 
(vhsqn)E (Xn, Mi) ||valli + llar llo 


: 定理 3.2 回 ”在 定理 2.1 和 定理 3.1 的 假设 条 件 下 , 问题 (13) 存在 唯一 的 解 (ws, pa) € (Xa, Ma) 
lurli < Tl, llpallo < 81 ( co 291110 + llfllo)- (18) 
4 ”误差 估计 
定义 如 下 的 Stokes 投 影 (Ra, Qh) : (X, M) > (Xn, Mn) 
Ba ((Ra(v, q), Qa (u, q)); (Vhs qn)) = B((o,q);i (vhs Gn), V (vhs gn) € (Xn, Mn). (19) 


314.15) 在 定理 3.1 的 假设 条 件 下 ， 对 任意 的 (wa) € (D(A), H1(Q) n M)，Stokes 投 
JÉ (Rn, Qn) 满足 


llv — Ra (v, allo + h( llv — Ra(v, alla + lla — Qal, allo ) < sh? (lole +llalli). (20) 
定理 4.2 在 定理 2.1 和 定理 3.2 的 假设 条 件 下 ， 成 立 如 下 估计 式 
lu —wuallo + h (llu — walla + lp — pallo) < Kh?. (21) 
证 明 由 (7) 减 (13) 得 


Bn ((en, Th); (vn, qh)) + bE + Eh, u, Un) a b(un, E+ eh, Un) =0, (22) 
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J E =u — Railu, p). s (22) 中 令 


(un, qh) = (en, m) = (Ra(u,p) — ux, Qalu, p) — ph), 


由 (3) 可 得 
vilen? + Glam) + blen ten) < p(B, €n)| + felun: E, e), 
这 里 
vlen? + B(en, nen) > v (1 — SÈ Jent, 
CE, u en)| + (uns E, en)| < of ul + unlla lea El < lleni: 
从 而 有 


llealla < Kh. 
类 似 地 ， 由 (13) 和 (7) 可 得 
Br ((en, Th); (Un, qh)) + b(u, E + eh, vn) 
+b(E + en, u, un) — b(E + ex, E + en, un) = 0. 
(vn, dn) = (ens mh)» FF 
vllenll? + G(m, mh) + blen, u, en) 
< |b(u, E,en)| + |b(E, u, en)| + |b(E, E,en)| + lb(en, E, en)|, 
再 由 (5) 和 (24)， 得 到 


lenll? < clb(u, B,en)| + |b(E,w, en)| + lb(E, E, en)| + blen, E, en)| 


A 


IA 


cllullallEllollealli + e( Ella + lezl )|Elhilles|i < pllenlla， 


即 
llenlla < sh’. 


ch (H1) 和 (30) 得 
lu — unllo < llenllo + Ello < rollenlla + «i?( lullz + lplli) < sh’, 
lu — walla < llex|li + El < sh? + sh( llullə + |lp|i ) < sh, 


利用 (6), (9), (15) 和 (16)， 得 


=I |b(u — up, u, vn) + b(un, u, va )| 

linalo < sup 上 
(wn,an)E(Xn,Mn) (llvalls 十 llaallo) 
< «(ul + llualls llu — walla < sllu — walls, 
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(23) 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


(30) 


(31) 


(32) 


(33) 
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4k 
4k 
ga 
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由 三 角 不 等 式 (30) 和 (33)， 得 
lp — pallo < lip — Qa (u, p)llo + lnnllo < sk( lulle + lplli) + sllu — ualli < rh. (34) 


H (20), (31), (32) 和 (34)， 即 得 (21) 成立 。 


5 ”数值 实验 


为 了 证 明定 理 4.2 中 收敛 率 的 正确 性 ， 本 节 考 虑 两 个 数值 算 例 。 稳 定 化 项 (12) 中 ， 局 部 压力 
投影 荆 的 选取 采用 了 Bochev 在 文献 [7] 中 的 技巧 。 对 所 有 的 数值 实验 ， 令 2 = [0,1] x [0,1], 
粘性 系数 > = 0.05。 给 定 符合 方程 (1)-(2) 的 真 解 如 下 


u= (ui(z,g)ua(zy))，p(z,y) = 10(z 一 0.5)， 


ui(x,y) = 10z2(1 — x)?y(1—y)(1—2y), uz(z,y) = —10z(1 — z)(1 — 2z)2(1 — y)?. 


f 由 方程 (1)-(2) 确定 。 
为 了 简单 起 见 ， 我 们 把 单位 正方 形 剖 分 成 N? 个 正方 形 元 或 2N? 个 三 角形 元 ，N 是 每 个 方 
向 上 等 分 的 数量 ， 步 长 h = 4. XHQ 一 Po TRMAR 7 TAA 
1 1 1 

h= g 16 — 5 56， 
X P, — Pa zú, 我们 选取 6 个 数值 

_ 1 1 1 1 

8 1624 48 
表 1 和 表 2 给 出 了 定常 Navier-Stokes 方 程 的 ?误差 和 五 1 误差 的 收敛 率 。 从 表 1 和 表 2 可 以 看 
出 ， 计 算出 的 收敛 率 十 分 逼近 第 4 部 分 理论 上 推导 的 数值 。 本 文 还 给 出 了 文献 占 提出 的 局 
部 稳定 化 方法 的 数值 实验 结果 ， 表 3 至 表 6 给 出 了 建立 在 三 角形 、 四 边 形 剖 分 基础 上 随 稳定 
化 参数 8 变化 的 数值 实验 结果 ;图 1 至 图 4 分 别 给 出 了 速度 相对 误差 ell 和 压力 相对 误 
Fe 号 给 + 随 稳 定 化 参数 变化 的 误差 曲线 。 显 然 ， 文 献 [5] 提出 的 局 部 稳定 化 有 限 元 算法 严重 依 
赖 于 稳定 化 参数 的 选取 。 但 本 文 应 用 的 局 部 稳定 化 有 限 元 算法 不 含有 任何 稳定 化 参数 。 总 之 ， 
稳定 化 有 限 元 算法 四 是 解决 定常 Navier-Stokes 方 程 的 一 种 有 效 的 数值 算法 。 


Ril: 压力 投影 稳定 化 方法 的 结果 : Q. — Po 


EE EE EENE 
E TES air 
0.884633 0.68003 / 0.162006 


1/24 


第 2 期 王 爱 文 等 ， 定 常 Navier-Stokes 方程 低 阶 混合 有 限 元 的 压力 投影 稳定 化 方法 255 


表 2: 压力 投影 稳定 化 方法 的 结果 : P, — Po 


u—u u—up l| H ak L 
lI Taela l hllz; rate lp 一 p 关 17 2 
2 lulla lpllz2 


/ 


0.835232 | / | 5682411 0.147155 
1.9008 0.158856 1.2938 0.0387494 
0.06370955 | 1.9347 0.111125 1.2421 0.0281026 
0.0412166 1.9516 | 0.0849661 1.2028 | 0.0220468 
0.0280234 2.1160 | 0.0686039 1.1732 0.0181435 


表 3: 稳定 化 方法 回 的 结果 : Qi 一 Po, B= 100 


alz Tola Tolz 
0.0704952 0.26942 0.125015 


0.0176969 1.9940 0.134964 0.9973 0.0625114 
0.0078708 1.9983 0.0900034 0.9992 0.0416761 
0.00442838 1.9992 0.0675092 0.9997 0.0312573 
0.00283447 1.9995 0.0540097 0.9998 0.0250055 
0.0019685 1.9997 0.0450091 0.9999 0.0208375 
0.00144629 1.9998 0.0385797 0.9999 0.0178603 


表征 ”稳定 化 方法 回 的 结果 : P, — Po, 8 = 100 
lu-urltz2 


lull z2 lulla llpll z2 
0.22867 / 0.483618 / 0.129906 


1/h 


[7 [sess [7 


AS: FRA! 的 结果 : Qi 一 Po, B = 0.0001 
luuria ETF ENA 
tales Mare | rate | e | 
1.05741 2.03577 0.255142 


1.01496 
1.01016 
1.00705 
1.00518 
1.00397 
1.00316 
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表 6: 稳定 化 方法 回 的 结果 :; P, — Po, B = 0.0001 


IPP Teller Piz 
0.999998 0.999998 0.628319 


-3.19701 
1.005806 
-2.41014 


oso | o | ossoa | -00099 [1.0880 _ 


-0.99916 


下 16*16mesheg 
-©- 24*24mesheş 
—+- 32*32meshes 


—— 16*16meshes 
—@— 24*24meshes 
一 | 一 32“32meshes 


图 1: 四 边 形 剖 份 回 的 速度 相对 误差 曲线 


—*— 16*16meshes 
~O- 24*24meshes 
~ 32*32meshes 


~+- 16*16meshes 
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图 3: 三 角形 剖 份 回 的 速度 相对 误差 曲线 ”图 和 ”三 角形 剖 份 回 的 压力 相对 误差 曲线 
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Stabilization of the Lowest-order Mixed Finite Elements Based on 
the Local Pressure Projection for Steady Navier-Stokes Equations 


WANG Ai-wen!, LI Jian?, XIE Dong-xiu' 


(1- Department of Mathematics, Beijing Information Science and Technology University, 
Beijing 100192; 2- Department of Mathematics, Baoji University of Arts and Sciences, Baoji 721007) 


Abstract: We present a stabilized finite element method based on the local pressure projection for 
the steady Navier-Stokes equations by the lowest order conforming finite element pairs (i.e. Qi — Po; 
P, — Po). In contrast to other stabilized methods, the new method is parameter free, it is not neces- 
sary to calculate higher order derivatives and edge-based data structures, and it is implemented at the 
element level with minimal cost. For the stabilized finite element solution, the optimal error estimates 
have been obtained. Finally, the numerical examples demonstrate the better stability and convergence 
of the stabilized finite element method for the Navier-Stokes equations. 

Keywords: Navier-Stokes equations; stabilized finite element; lowest order pairs; local pressure pro- 


jection 
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